Large-Scale Structure with Gravitational Waves II: Shear 
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The B- (curl-)mode of the correlation of galaxy ellipticities (shear) can be used to detect a stochas- 
tic gravitational wave background, such as that predicted by inflation. In this paper, we derive the 
tensor mode contributions to shear from both gravitational lensing and intrinsic alignments, using 
the gauge- invariant, full-sky results of Schmidt and Jeong p|. We find that the intrinsic alignment 
contribution, calculated using the linear alignment model, is larger than the lensing contribution 
by an order of magnitude or more, if the alignment strength for tensor modes is of the same order 
as for scalar modes. This contribution also extends to higher multipoles. These results make the 
prospects for probing tensor modes using galaxy surveys less pessimistic than previously thought, 
though still very challenging. 



I. INTRODUCTION 

A stochastic gravitational wave (GW) background is 
one of the key testable predictions of inflation. Thus, 
much theoretical and experimental effort is devoted to 
searching for this gravitational wave background. One 
necessary ingredient of any method designed to search for 
a GW background is the ability to cleanly separate the 
GW contribution from scalar perturbations. Most com- 
monly, this is done by considering spin-2 quantities on 
the sky, such as the anisotropy in the polarization of the 
cosmic microwave background (CMB) radiation 
the anisotropy of the 2-point correlation function 
or the ellipticity of galaxy images Q. The polarization 
of the cosmic microwave background is commonly con- 
sidered as the most promising probe. However, the 21cm 
emission from the dark ages has recently been shown to in 
principle offer even more discovery potential 0, Q . On 
the other hand, previous authors have concluded that 
weak lensing shear will most likely not be a competitive 
probe of primordial GW @, EH- Nevertheless, given 
the scientific impact, it will be crucial to confirm a pos- 
sible detection of a GW background in the CMB with 
independent methods, such as shear. 

The goal of this paper, and its companion pj, is to 
systematically and rigorously derive the GW effects on 
large-scale structure observables. While we restrict our- 
selves to a linear treatment in the tensor perturbations, 
we strive to keep the results as general as possible other- 
wise. 

This paper deals with shear, i.e. the correlations of 
galaxy ellipticities. The underlying assumption in inter- 
preting shear correlations is that in the absence of per- 
turbations, galaxy ellipticities are uncorrelated on large 
scales. However, in the relativistic context, this raises the 
question of the frame of reference in which galaxy ellip- 
ticities are actually uncorrelated. In Schmidt and Jeong 
we have derived an expression for the shear based 
strictly on observable quantities ("standard rulers"). As 
shown there, that expression is equivalent to the state- 



ment that galaxies are isotropically oriented in their local 
inertial frame, described by the Fermi normal coordinates 
along a given galaxy's geodesic. 

Consider a region of spatial extent R, much larger than 
the size of individual galaxies but smaller than the typ- 
ical wavelength of the perturbations we aim to measure 
through shear. The center of mass of this region moves on 
a time-like geodesic. We can construct a coordinate sys- 
tem where the center of mass is at rest at the origin, and 
the time coordinate tp corresponds to the proper time 
of this geodesic. In other words, the unit vector defin- 
ing the time coordinate is equal to the tangent vector to 
this geodesic. The spatial coordinate lines are chosen to 
be space-like geodesies ("straight lines") orthogonal to 
this time direction, and whose unit vectors are parallel- 
transported along the observer's geodesic. These are the 
so-called Fermi normal coordinates (FNC) pj|, [l3j]. It 
is straightforward to construct them perturbatively for 
a given metric and time-like geodesic (see App. [A}. The 
most important property of Fermi normal coordinates is 
that the metric is Minkowski at the spatial origin at all 
times tp, with corrections away from the geodesic go- 
ing as x F x F jR 2 c , where x F denote the spatial Fermi co- 
ordinates and R c is the typical curvature scale of the 
space-time (R c ~ H^ 1 for an unperturbed FRW space- 
time). Thus, the Fermi normal coordinates are the frame 
in which a local observer in a weak gravitational field 
would describe her experiments. Neglecting the correc- 
tions oc x F , there is no preferred direction in these coor- 
dinates along which galaxies could align, and they thus 
have to be oriented isotropically. The contribution to the 
shear from the transformation to the local Fermi coordi- 
nates was first introduced by Dodclson et al. Q, who 
showed that this term, which they call "metric shear", 
has a significant impact on the observational signatures 
of tensor modes. Furthermore, as shown in p|, this term 
is crucial to ensure that the expression for the observ- 
able shear does not receive contributions from constant 
and pure-gradient metric perturbations, as required by 
coordinate-invariance. 

Here, we evaluate this expression for the shear in an 
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FRW space-time with tensor modes, which yields the 
shear contribution by "projection effects" induced by a 
GW background. As discussed in 0, EH , the tensor con- 
tributions typically peak at the source location rather 
than at lower redshifts, and there is no enhancement of 
the contribution by transverse modes. Both of these facts 
lead to qualitative differences from the scalar case. 

As discussed, the space-time in Fermi coordinates 
around a given galaxy is not perfectly flat however, and 
the corrections to the metric oc x l F x F /R 2 provide pre- 
ferred directions along which galaxies can align. For 
non-relativistic motions (which generally applies to large- 
scale structure), the relevant contribution to the metric 
in the Fermi frame corresponds to a tidal field, which can 
align galaxies and thus contribute to the observed shear 
correlation. Here we derive the contribution to the tidal 
field by tensor perturbations, and for the first time calcu- 
late the intrinsic alignment contribution of tensor modes 
in the linear alignment model. This prescription has been 
shown to agree well with observations on large scales in 
the scalar case. Note that this approach is only appli- 
cable for gravitational waves with periods much longer 
than the dynamical time of galaxies. 

The key advantage of the shear applies to both lens- 
ing and linear intrinsic alignment contributions: linear 
scalar perturbations contribute only to the parity-even 
i?-(gradient-)mode component at linear order, while ten- 
sor perturbations also contribute to the B-(curl-)mode 
(both through lensing and intrinsic alignments). Impor- 
tantly, scalar perturbations do contribute to i?-modes at 
second and higher order, a point which we will discuss in 
detail in § IVTAl 

The outline of the paper is a follows: we introduce our 
notation and conventions in § |TT] (they are the same as 
in the companion paper Q ) • § Mil presents the deriva- 
tion of the lensing (projection) contribution, while § IIVI 
discusses the intrinsic alignment contribution from tensor 
modes. §IVl gives the expressions for the E— and -B-mode 
power spectrum of the shear, including the connection to 
previous results. § I VII presents the results. We conclude 
in § IVIII The appendix contains details on Fermi normal 
coordinates, the derivation of angular power spectra, and 
the connection to convergence and rotation. 



II. PRELIMINARIES 

We begin by introducing our convention for the metric 
and tensor perturbations and some notation; it is iden- 
tical to that used in 0. For simplicity, we restrict our- 
selves to a spatially flat FRW background, and consider 
only tensor (gravitational wave) modes in the main part 
of the paper. The perturbed metric is given by 

ds 2 = a 2 (i]) [-dry 2 + (<% + hy) dx l dx 3 ] , (1) 

where hij is a metric perturbation which is assumed to 
be transverse and traccless: 



We then decompose hij into Fourier modes of two polar- 
ization states, 

MM) =e+(k)/ l +(M) + e*(k)h*(k, V ), (3) 



where ef.-(k), 



x , are transverse (with respect to 



k) and traccless polarization tensors normalized through 



25 s 



Note that h s = e\ J hij/2. Wc assume 



both polarizations to be independent and to have equal 
power spectra: 

(h s (k, V )h s , (k', 77')) = (2n) 3 S D (k - k')S ss , ip T (fc, T], i!'). 

(4) 

Here, 77 denotes conformal time, and the unequal-time 
power spectrum is given by 



P T (k, 77, r/) = T T (k, r,)T T (k, 7/)P T0 (fc), 



(5) 



where T^(/c,?7) is the tensor transfer function, and the 
primordial tensor power spectrum is specified through 
an amplitude A T and an index tit via 



P T0 (k) =27T 2 fc" 



(6) 



Following WMAP convention (l5| . we choose ko = 
0.002 Mpc -1 as pivot scale. Throughout, we will assume 
a scalar-to-tensor ratio of r = 0.2 at ko (consistent with 
the 95% confidence level WMAP bound), which together 
with our fiducial cosmology determines A T . The tensor 
index is chosen to follow the inflationary consistency re- 
lation, tit = —r/8 = —0.0025. For the expansion history, 
we assume a flat ACDM cosmology with h ~ 0.72 and 
n m = 0.28. Contributions from scalar perturbations are 
evaluated using a spectral index of n s = 0.958 and power 
spectrum normalization at z = of as = 0.8. 
From Eq. ([3]) and Eq. (TJ}, we easily obtain 



(h^k, r))h kl (k', 77')) = C2n) 3 5 D (k - k') 



(7) 



(k)e+(k) + e>Kk)e*(k) 



x ^-Pr(fc, 77,77') 

{hij (k, 7,)^'(k', 7/)) = (27r) 3 S D (k - k')P T (k, 77, 77'). 

Long after recombination, the transverse anisotropic 
stress which sources gravitational waves becomes neg- 
ligible, and the tensor modes propagate as free waves. 
During matter-domination, the tensor transfer function 
then simply becomes 



Tt(M) 



krj 



(8) 



which however is still valid to a high degree of accuracy 
during the current epoch of acceleration. Wc will use 
Eq. ([5]) throughout. We also define 



= (h lk y 



(2) 



pij = S ij _ „',v 



(9) 
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as the projection operator perpendicular to the line of 
sight. 

As a traceless 2-tensor on the sphere, the shear can 
be decomposed into spin±2 functions ±27 = 7i ± 272 (in 
analogy to the combination of Stokes parameters Q ±iU 
for radiation) as 

lij = 2l'm\m 3 j r + _27m'_m 3 _. (10) 

Here, we have defined the unit vectors of the circularly 
polarized basis, m± = (ee^fie^)/^ (see App. A in [l[). 
If we choose a coordinate system where h l is along the 
z-axis and e l g along the x-axis, we have 

(7i 72 \ 
72 -71 . (11) 
/ 

This decomposition is particularly useful for deriving 
multipole coefficients and angular power spectra. In 
particular, we can define the multipole moments of the 
parity-even .©-modes and parity-odd S-modes through 

a Jm = 5 ( a 7m + a Tn) 

«7™ = - - a Z) ■ (12) 

We also define 

X± = m^X t 

E± = rri^m^Eij (13) 
for any vector Xj and tensor E\j. 

III. LENSING EFFECTS 

We begin with the definition of shear in the relativistic 
context as derived in [l[, for a synchronous comoving 
metric as in Eq. ([T]): 

±27 = m^m^Aij 

= - -h± - m l T m : > T d ±i Ax_ Lj . (14) 

In other words, the shear is the traceless part of the sym- 
metric 2x2 matrix Aij which describes the transverse 
distortion of transverse standard rulers. Ax l ± is the dis- 
placement perpendicular to the line of sight of the ob- 
served position from the true position of the source, in 
terms of the global comoving coordinates. 

As shown in [l[, Eq. (fT4| is explicitly given by 

{±2l){ti) = -\h± -\h ± (15) 
- f X dxl^^-irn^m^didjhki^fi 1 

Jo y 1 x 

+ h l m%m%dih u - -U± I. 



The second term on the r.h.s. can be understood as com- 
ing from the transformation from global coordinates to 
the local Fermi normal coordinates [l[; in the follow- 
ing we will refer to this as the "FNC term". It is im- 
mediately clear that a constant metric perturbation hij 
(which corresponds to a pure gauge mode) does not con- 
tribute to 71 ± ij2- The same is true for a pure gradient 
hij = BijkX k . In App. C of [l| we have applied several 
test cases to Eq. (|15l) . including a Bianch I cosmology 
where all terms contribute non-trivially. 

IV. INTRINSIC EFFECTS 

Eq. (fl~5|) is derived assuming we have a perfect stan- 
dard ruler, in the sense that the intrinsic physical size 
of the ruler is uncorrclated with the perturbations hij. 
In the case of weak lensing shear surveys, the "standard 
ruler" is the fact that galaxies are randomly oriented, i.e. 
their apparent size measured along two different fixed 
directions is on average equal. However, we know that 
in reality galaxy orientations are not truly random, but 
there is some alignment with large-scale tidal fields. 

In order to determine the effective tidal field experi- 
enced by galaxies in a Universe with propagating ten- 
sor modes, we derive the corrections cx x 2 F to the metric 
g F v in the Fermi normal coordinate frame. In particular, 
since we are concerned with non-relativistic motions, we 
are mostly interested in the time-time part of the metric 
g F . The detailed derivation for a space-time described 
by Eq. (TT|) is presented in App. [A] The result for g F is 

g F (x Fl t F ) = -\ + (H + H 2 )r 2 F -2* F {* F ,t F ) (16) 

* F (x F ,£ F ) = - -[hij(0,tF)+2H(t F )hij(0,t F )]x F x F , 

(17) 

where r F = SijX F x F , and dots indicate derivatives with 
respect to time t (equivalent to t F at this order). The 
terms cx r F in Eq. (|16[) arc the usual Hubble drag which 
provide an effective repulsive force. The leading effect 
of large-scale cosmological perturbations on the region 
considered is to add an effective tidal field fy F , which 
depends on the time derivatives of the metric perturba- 
tion hij. Note that for a traceless tensor perturbation, 
y F is indeed a purely tidal field (i.e. V 2 * F = 0). The 
fact that the amplitude of the tidal field is given by the 
time derivatives of the tensor modes implies that only 
modes that have entered the horizon contribute; specifi- 
cally, the contribution for a given Fourier mode scales as 
k 2 for k — > 0, as can be seen from Eq. §8§ and Eq. (|2"0")) 
below. 

In the absence of perturbations (hij = 0), Eq. p^|) 
is isotropic. Thus, in the Fermi frame, there is no pre- 
ferred direction along which galaxies forming in this re- 
gion could align, and their orientations are truly random 
in this case. As shown in [![, the expression for the shear 
Eq. (|14[) contains precisely this statement (of course this 
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holds not only for galaxy orientations, but any standard 
ruler). 

In the presence of perturbations, the tidal field 
provides a preferred direction along which galaxies can 
align. The fact that large-scale tidal fields tend to align 
galaxies (intrinsic alignment, IA) is well established both 
theoretically and observationally for scalar perturbations 
(e.g., [la-yjl)- In order to make progress, we will adopt 
the linear alignment (LA) model, which has recently been 
shown to be consistent with observations on large scales 
(> 10Mpc/M [ll,[2(|. In this model (following the no- 
tation of [I!]!]), the tidal tensor iy at the location of the 
galaxy, defined through 

tij = (didj - ^<%V 2 ^) *, (18) 

contributes to the traceless part of the observed distor- 
tion matrix Aij of the galaxy image through 

= - \c lPcM H- 2 V lk V 3 it kl {z p ). (19) 

That is, A\f is proportional to the projection of the tidal 
tensor onto the sky plane. Here, p CI o = 3Hq/(8-kG) is the 
critical density today. The constant of proportionality 
C\ determines the magnitude of alignment, while z p , the 
redshift at which the tidal field is evaluated, is another 
parameter of the model. Observationally, C\p cr Q ^0.1 
for galaxies at redshifts less than 1, when choosing z p 
to be equal to the source redshift. A positive C\, to- 
gether with the overall sign in Eq. (fT9|) . corresponds to a 
galaxy's major axis aligning with the smallest eigenvec- 
tor of Uj] physically, an initially spherical perturbation 
will tend to collapse last in the direction of the smallest 
potential curvature, leading to a preferential alignment of 
the major axis with this direction. Note that in Eq. (|T9"1) . 
Cip cr o multiplies the tidal field in physical rather than 
comoving units. 

While this mechanism is expected to be qualitatively 
the same for tensor modes as for scalar modes, there is no 



reason for the amplitude C\ to be the same in both cases. 
In particular, linear tidal fields sourced by scalar per- 
turbations are constant during matter domination, while 
tensor perturbations decay and oscillate. However, we 
will see that the tensor modes relevant for the IA contri- 
bution arc long-wavelength and should not have a qual- 
itatively different impact on the formation of galaxies 
than scalar tidal fields. We will return to this issue in 
§ IVII In the following, we will assume that z p is equal 
to the source redshift. Generally, evaluating the tidal 
field at z p > z leads to larger effects so that this is a 
conservative assumption. 

Using Eq. (jTTJ) and Eq. (| 14[) . it is then straightforward 
to evaluate the contribution of tensor modes to the shear 

±27 : 

( ±27 ) IA (n)=m^4^ A 

= ^C lPci . H- 2 (h± + 2Hh± S ) 

= \Ci Pcl0 H^ 2 a- 2 (hi + aHh' ± ) , (20) 

where primes indicate derivatives with respect to 77, and 
we have used dt = adr\. 



V. OBSERVED SHEAR STATISTICS 

We now use the results derived in the previous two 
sections to calculate the angular power spectrum of the 
observed shear induced by tensor modes. We briefly out- 
line the steps of the derivation, which follows the general 
prescription described in App. Al of [l|, with details rel- 
egated to App. [Bj 

In the first step, we consider a single tensor mode of 
wavevector k which we assume oriented along the z-axis. 
Including the intrinsic alignment contribution (§ II V[) . we 
obtain for the contribution to the shear 



( ±27 )(k,n)= J2 \-\ 
p=-i,i I 



- ^C lPcr0 H^ 2 d- 2 {d? +aHdfi }j h p (k,f,)e lk& < 



-(lTPAt) e 



2 i2p0 



i_^| fc2(1 _ + _ 2 |) _ M a )(1 T w) + J_ (1 T mf 



where h±\ = (h + =p ih x )/2 are the circular polarization states of tensor modes, /i = cos#, and 9, <p denote the polar 
and azimuthal angles of the line of sight unit vector n. We then apply the spin- lowering/raising operator twice to ±27 
in order to obtain a rotationally invariant (spin-0) quantity. The spherical harmonic coefficients of this scalar quantity 
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are directly related to those of the spin±2 shear components (see App. A in [lj). We can then separate a( m into a 
real part a] E , which is parity-even and thus transforms in the same way as a tensor derived from a scalar function /, 
and an imaginary part io,7^, which acquires an additional minus sign under parity. This E/B- mode decomposition 
is useful because any symmetric tensor 7^ derived from a scalar function is parity-even and thus does not source any 
.B-modes (as shown explicitly for the shear in Further, any perturbations generated by parity-conserving physics 
do no induce an E — B cross-correlation. 

Finally, using relations derived in App. Al of we obtain the angular power spectra of E- and B-modes of the 
shear induced by tensor modes (App. IB 2j) : 



c xx {i) 



F? E (k) 



F? B (k) = 7 



i- J k 2 dk P T0 (k)\F? X (k)\ 2 , X = E,B 



(21) 



Tt(Mo) RcQi(x 



ReQ 2 (ir) + ^ReQ 3 (a:) 
X 



x dx 
X 

T T (k, 770) f ImQi(j;)^- 



,r=0 

3l{x) 



1 - ^C lPcl0 Ha 2 a- 2 [d 2 + aHdrM T T {k, ry)Re &(a 



T T (k,r] - x) 



x=0 



1 — gCip cr oi? 



" 2 a" 2 



{dl + aHdfj}) T T (fc,77)ImQi(x) 



X x 



Here Qi(x) are derivative operators whose action on 
ji(x)/x 2 is given explicitly in Eq. (IB16[) : in particular 
lmQ 3 (x)[ Jl (x)/x 2 ]=0. 



A. Relation to convergence and rotation 

The shear is most commonly written in terms of an- 
gular derivatives of the deflection angle A8\ i.e. 7^ is 
defined as the trace- free part of dA6j/d0i. In our nota- 
tion, A8 l = Ax l ± /x, and dA8 :j /dd l = d ±i Ax ±j . There 
are two degrees of freedom in Ax l ^_ , and hence only two 
independent components of d±iAx±j. If we define the 
convergence k and rotation u> through 



8 ±i Ax i 1 



u = ~ 2 £ ijkn l d 3 L Ax'i, 



(22) 



then the i^-mode of the shear defined as trace-free part 
of d±iAx±j is directly related to k, while the J5-mode is 
related tow. As shown by Stebbins [2l|, on the full sky 



,EE m _ + 2)0-1) 
"> 1(1 + 1) 



C ft (Z) 



(23) 



However, neither of the deflection angle Ax* ± , the con- 
vergence k, or the shear defined as trace-free part of 
d±iAx±j are observable (the rotation is only observable 
if there is an intrinsic preferred direction in the source 
plane). Instead, the observable shear is given by Eq. (fT5)) 
which includes the FNC term. Nevertheless, the relations 



Eq. (|23[) arc useful as an analytical and numerical cross- 
check of Eq. (f2"Tj) (without the IA contribution), and this 
cross-check is presented in App. O 

In the previous work of Q , the rotation was used as a 
proxy for shear _B-modcs through Eq. ((23]) . This result 
thus does not include the FNC term. However, they also 
consider the "metric shear" , which is the contribution 
to lu that corresponds to the shear contribution by the 
FNC term. When including the metric shear in w, we 
indeed recover the B-modes of the shear including the 
FNC term through the relation Eq. fl55J} (see App. EH). 
Thus, our results for the lensing-induced shear £?-modcs 
from tensor modes (i.e., neglecting the intrinsic align- 
ment contribution) agree with those of 0, modulo the 
factor (l + 2)(l-l) 



VI. RESULTS 

We begin by investigating the separate terms con- 
tributing to the shear E- and -B-mode power spectra in 
Eq. (|2ip . focusing on the lensing contributions without 
intrinsic alignment first. Fig. Q] shows the -E-mode an- 
gular power spectra obtained by separately considering 
only the terms cx Q2, oc Q3, and the FNC and observer 
terms cx Q\. Here, we have multiplied the power spec- 
tra by I 6 as they are very steeply falling with I. We see 
that for I < 10 there is a very significant cancellation (by 
three orders of magnitude) between the different terms. 
In fact, the magnitude of the individual terms depends 
on the lower limit used in the integration over k, and di- 
verges logarithmically as k m \n — > for small I. For the 
B-modes, the cancellation is not as important, though 
still significant (Fig. [5] in this case there is no Q3 term). 
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FIG. 1: Lensing (projection) contributions to the observed an- 
gular power spectrum of the E-mode component of the shear 
from tensor modes, separated into terms oc Q± (observer and 
FNC terms), oc Qi and oc Q3 respectively. Note that power 
spectra are multiplied by I 6 . We have assumed a sharp source 
redshift of z — 2. 



This result is in agreement with the findings of [8[ : when 
including the FNC contribution (metric shear), the am- 
plitude decreases significantly. This is not surprising: if 
we drop the FNC term in Eq. (|2ip. tensor perturbations 
with k — > contribute to the shear at low I, and since 
the tensor power spectrum is sharply falling with k, these 
contributions are large. Clearly, these contributions are 
unphysical however, and the FNC term must be included. 

Fig. [3] shows the total lensing contribution to the E- 
mode power spectrum of the shear, the intrinsic align- 
ment contribution, as well as the sum of the two, while 
Fig. U shows the same for the (more interesting) B- 
modes. Here, we assumed a Gaussian redshift distribu- 
tion dN/dz oc exp(-(5 - z) 2 /2Az 2 ) with Az = 0.03(1 + 
z). Further, we have adopted a value of Cip rT n = 0.12 
as measured in the Sloan Digital Sky Survey [ljj. This 
coefficient will depend on the specific galaxy sample con- 
sidered, in particular on the redshift. Here we extrapolate 
the value of Cip cr o measured for galaxies at z « 0.3 — 0.5 
to galaxies at z = 2, assuming no evolution. Thus our 
results should only be seen as a rough estimate of the 
magnitude of this effect (note however that we assume 
a constant alignment strength with respect to the phys- 
ical, not comoving tidal field). Even with this caveat 
in mind however, it is clear that the intrinsic alignment 
contribution is far larger than the lensing contribution. 
This is in contrast to the scalar case, where for source 
galaxies at cosmological redshifts the lensing signal is 



FIG. 2: Same as Fig. [2l but for the B-mode component. 



significantly larger than the intrinsic alignment contri- 
bution. The underlying reason is that the projected con- 
tributions from lensing arc relatively suppressed in the 
tensor case. While scalar perturbations with transverse 
wavevector deflect light coherently along the past line 
cone to the source, tensor perturbations propagate and 
decay, such that no such coherent deflection occurs even 
for transverse wavevectors [3]. The result is that the 
lensing contributions are mostly localized at the source 
for tensor modes, and down- weighted by the lensing ker- 
nel (oc X-x). 

Apart from Ci/9 cr o, the linear alignment model has an- 
other free parameter in the redshift z p at which the tidal 
field is evaluated. By default, we choose z p = z. How- 
ever, choosing z p to correspond to a time 5 x 10 s years 
before observation (z p sa z + 0.04 for 2 = 2), which cor- 
responds to several dynamical times for typical galax- 
ies, only yields a slight increase in the power spectrum 
contribution by ~3%. Varying z p thus does not have 
a significant impact on the intrinsic alignment contribu- 
tion. On the other hand, this mild dependence on z p 
indicates that the bulk of the IA contribution induced 
by tensor modes is due to slowly varying tidal fields, i.e. 
tensor modes with k/H ~ 1, rather than rapidly oscillat- 
ing modes with k/H 3> 1 (this is confirmed by numeri- 
cal inspection of the intrinsic alignment contribution to 
F^' B ). Such tidal fields, which vary on a Hubble time, 
are not expected to have a qualitatively different effect 
on the formation of galaxies and halos than scalar tidal 
fields, given that the relevant time scale is the dynam- 
ical time of the collapsing dark matter halo. We thus 
expect that the value of Cip C ro relevant for the I A con- 
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FIG. 3: Angular power spectrum of the observed E- 
component of the shear from lensing and intrinsic alignment 
affects, as well as the total power spectrum. We assumed 
Cipci-o = 0.12 (following the results of and a Gaussian 

distribution of source redshifts centered at z = 2 with RMS 
width of Az = 0.03(1 + z). 



tributions to shear from tensor modes will not be very 
different from that for scalar tidal fields. However, one 
would expect Ci/9 cr o to be generically scale-dependent for 
tensor modes, decaying from its low-fc limit to smaller 
values as 1/fc approaches the scale of halos and galaxies 
(fc > 0.3/i/Mpc). 

Fig. [5] shows the redshift evolution of the -B-modes of 
the shear. As expected, larger source redshifts yield sig- 
nificantly larger signals, due to the decay of the tensor 
modes and since at higher redshifts, larger scales are be- 
ing probed at a given I. However, we also see that the 
intrinsic alignment contribution evolves even faster with 
source redshift (note that here we have assumed the same 
value for Cip cr o at all redshifts). This can be traced back 
to the factor of a~ 2 in the IA contribution [Eq. (j^Tj) ] . 
which is due to the transformation from conformal time 
derivatives to physical time derivatives. It is also interest- 
ing to consider the dependence of the signal on the width 
of the source galaxy redshift distribution. This is illus- 
trated in Fig. [6l The lensing contributions are largely 
independent of Az for the range of multipolcs relevant 
here. On the other hand, the IA contribution is notice- 
ably increased for sharp source redshifts at I > 10, a 
consequence of the fact that this contribution is not pro- 
jected along the line of sight but evaluated at the source. 
Thus, unlike the lensing contribution, the IA contribution 
is essentially a three-dimensional field. Note also that in 



FIG. 4: Same as Fig. O but for B-modes. 



this case 1(1 + l)C^ B (l) w const, i.e. there is roughly 
equal power per decade in multipole for the IA contribu- 
tion. However, following our discussion above, we expect 
the approximation of a scale- independent alignment coef- 
ficient to break down once the wavelength of contributing 
tensor modes approaches the scale of halos, roughly at I 
greater than a few hundred. 



A. B-modes from scalar perturbations 

In addition to being sourced by primordial tensor 
modes, shear B-modes are also produced by second-order 
corrections to lensing by scalar perturbations (here we ne- 
glect second order corrections to the scalar intrinsic align- 
ment contributions). These come from three sources. 

The first source is tensor modes generated by non- 
linear gravitational instability (e.g., [13, H3l)- The shear 
-B-mode power spectrum induced by these tensor modes 
was found to be roughly scale-invariant and at the level 
of 1(1 + 1)C bb (1)/2tt ps 1CT 14 Q3 which is much smaller 
than the total primordial GW signal presented here for 
r = 0.2. However, the calculation of [10( did not include 
the intrinsic alignment effect which will also increase the 
B-mode signal of scalar-induced tensor modes. Due to 
the different scale-dependence and redshift evolution of 
the latter, the boost will likely be somewhat smaller than 
that for the primordial tensor modes. We leave this for 
future work. 

The second, more significant scalar source for i?-modes 
is from second-order corrections in the geodesic equation 
(beyond-Born correction and lens-lens coupling) (24 . HBJ . 
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Wc have evaluated these according to 



cBB ^(l + 2)(l-l) 



1(1 + 1) 



I 7 



(24) 



where Cf u is given in Eq. (51) of [25|. Note that the 
expression given there uses the Limber and flat-sky ap- 
proximations and is thus not expected to be accurate 
for I < 10. Wc compare this contribution with the to- 
tal contribution from primordial tensor modes (intrinsic 
alignment and lensing) in Fig. [Jj At z = 2, the scalar 
contributions become larger than the tensor mode signal 
at I > 6. For higher redshifts, they only dominate at 
higher I. 

Finally, a third contribution comes from the fact that 
observationally we measure the reduced shear g = 7/ (1 — 
k) [2(||27|. Furthermore, selection effects ( "lensing bias" ) 
lead to a similar second order correction [28| . At leading 
order, both contributions can be summarized by writing 
the observed shear tensor as 



7 ° bs (n)= 7y (n) + (l + g Hn) 7li (n), 



(25) 



where the parameter q parametrizes the lensing bias con- 
tribution [2g. We can evaluate the scalar S-mode con- 
tribution from both of these effects using Eq. (21) in [28j . 
Note that this equation was derived in the flat-sky limit 
and hence will also not be accurate at I < 10. This con- 
tribution is also shown in Fig. [7J (assuming a lensing bias 
coefficient of q — 1). This contribution is even larger 
than that from the second order Born correction, and 
in fact dominates over the primordial GW contribution 
at z — 2 for all I but 1 = 2. In principle one could re- 
duce this contribution significantly by selecting a source 
galaxy sample with g w —1, although whether this is 
feasible in practice would need to be investigated. The 
fact that the reduced shear and lensing bias contribu- 
tions produce the dominant scalar contribution to shear 
B-modes is an interesting result in itself and has not been 
pointed out before. 



VII. DISCUSSION 

In this paper, we have studied the shear induced by a 
primordial GW background. In addition to the projec- 
tion (lensing) effects, for which we use a gauge-invariant 
expression, we derive for the first time the contribution 
due to intrinsic alignment (IA) of galaxies through the 
effective tidal field induced by tensor modes. We have 
found that this contribution is typically much larger than 
that from lensing. While surprising initially, this is due to 
the qualitatively different properties of lensing by tensor 
modes as compared to scalar modes. 

The IA contribution depends on a coefficient which 
can be observationally determined for scalar perturba- 
tions. In general, this does not have to be the same for 
tensor modes, since tensor modes evolve while the scalar 
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FIG. 5: Dependence of the lensing and intrinsic alignment 
contributions to the B-mode shear power spectrum on the 
source redshift z. We have assumed a Gaussian redshift dis- 
tribution centered at z = 5, 2, 1 (from top to bottom), and 
RMS width Az = 0.03(1 + z). The black dotted line near 
the top of the figure shows the la error on the shear power 
spectrum per multipole induced by shape noise [Eq. (|27[) ]. for 
a survey with n = 100 arcmin -2 , a e = 0.3, and / s k y = 0.5. 



tidal field is constant on large scales (during matter dom- 
ination). On the other hand, the bulk of the tensor con- 
tributions is from horizon-scale modes, which evolve on 
the Hubble time scale. Compared to the dynamical time 
of galaxies and halos, this is a very slow evolution, and 
we expect the alignment coefficient to only be mildly af- 
fected by this. On the other hand, the results shown in 
Fig. [3] through [5] depend on the alignment strength at 
high redshifts, which is currently poorly known observa- 
tionally. 

The IA contribution also decays much more slowly to- 
wards high I than the lensing contribution (especially for 
narrow source redshift distributions). In principle, this 
could allow one to access smaller-scale tensor modes than 
those probed by the cosmic microwave background. Fur- 
ther, since the IA contribution is not projected along the 
light cone, this effect in principle allows us to measure the 
entire three-dimensional field of tensor perturbations. In 
the case of perfect redshift measurements, the number of 
tensor modes measurable up to a maximum scale would 
thus scale as rather than 

On the other hand, even with this increased signal, the 
requirements for a detection of a stochastic GW back- 
ground using galaxy ellipticities are still extremely chal- 
lenging. Purely in terms of statistical power, the intrinsic 
ellipticities of galaxies add noise ( "shape noise" ) to the 
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FIG. 6: Dependence of the lensing and intrinsic alignment 
contributions to the B-mode shear power spectrum on the 
width of the source redshift distribution Az. z — 2 for all 
curves. 



shear E- and _B-modc power spectra. The corresponding 
la error on the power spectrum per multipole is given by 



(0 



1 



y/(2l + l)/ sky 2n 



(26) 



where X = E,B, / s k y is the fraction of sky covered by 
the survey, a e is the RMS intrinsic cllipticity of galaxies, 
and n is the number of source galaxies per steradian. In 
numbers, this yields 



l l±^-AC^ x (l)^1.6x 10- 11 , 
2tt 7 w V 2 



7 \ 3/2 

l\ ,-1/2 ( ^ 
/sky 1 n o 



0.3 



100 arcmin 



(27) 



This prediction, for / s k y = 0.5, is shown as black dotted 
line in Fig. [5j Clearly, one would need to go to source 
rcdshifts z > 2, at very high source densities, to detect 
a GW contribution at the level of r = 0.2 (unless the 
alignment strength at z ~ 2 is significantly larger than 
that measured at low redshifts). 



Given the smallness of the signal and the possible con- 
tamination by scalar contributions (S I VI Ap . quantitative 
constraints on an inflationary GW background will thus 
be very challenging. Nevertheless, shear S-modes remain 
one of only a handful of probes in cosmology that can be 
used to search for such a background. 



10" 



CQ 



io- 



10" 



CQ 

u 



10- 



10" 



io- 




Tensor 

Scalar (Born) 

Scalar (red. shear) : 

I ,1,1,1 I , i I , I n 

4 6 8 10 20 40 6080 
I 



FIG. 7: Shear B-modes from tensor perturbations (black 
solid, including both lensing and intrinsic alignment) and sec- 
ond order scalar contributions from corrections to the Born 
approximation (magenta dashed) and reduced shear and lens- 
ing bias (assuming q = 1, red dash-dotted). Results are for 
redshifts z = 5, 2, 1 (from top to bottom). Note that the 
scalar contributions have been calculated using the Limber 
and flat-sky approximations, respectively, which are not ac- 
curate at I < 10 (dotted lines). 
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Appendix A: Fermi normal coordinates 



In this appendix we review the basic concept of Fermi normal coordinates and its application to a flat FRW metric 
with perturbation hij of the spatial components (Eq. (UJ , but without imposing the transverse or tracefree conditions 
on hij). General covariance allows us to choose coordinates such that, at any given space-time point P, the metric is 
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Minkowski and the Christoffel symbols vanish [29[; i.e. 



9p-v — Vtiv 

9fJ.f,a 0. 



(Al) 



Riemann normal coordinates realize such coordinates by using a set of four geodesies (one time-like and three space- 
like) starting from the fixed point P. Fermi normal coordinates (FNC) are a specific extension of Riemann coordinates 
such that Eq. (|A1|) holds for every point along a fixed time-like geodesic. Specifically, we single out a time-like geodesic 
("central geodesic") that passes through P, as the worldline of the observer around which the FNC are constructed. 
Given three space-like tangent vectors at P which are orthogonal to the tangent vector of the central geodesic at 
P, the tangent vectors at all other points along the central geodesic are defined through parallel transport. For all 
points along the central geodesic, we construct Riemann coordinates using these tangent vectors. Then, the condition 
Eq. (|A1[) is satisfied at all points along the central geodesic. Given a central geodesic, FNC are uniquely defined up 
to three Eulcr angles. The significance of FNC is that they are the natural coordinates in which an observer moving 
along the central geodesic would describe local experiments. With the conditions Eq. (|A1[) satisfied, the departure 
from Minkowski of the metric in FNC appears at quadratic order in the spatial Fermi coordinate x l F . 

The Fermi normal coordinates can be explicitly constructed as follows (see [HI, [3(3]). The time coordinate is chosen 
to coincide with the proper time along the central geodesic. Let x F (P) = tp denote its value at point P. We can 
construct the spatial slicing of FNC, i.e. the x° F = const hypersurface, as comprising all space-time points in the 
neighborhood of the central geodesic that are reached by a congruence of spatial geodesies whose tangent vector at P 
are orthogonal to the tangent vector of the central geodesic (eo) P . Let Q denote a point on this hypersurface, so that 
x a F (Q) = tp. Further let x^(X) be the unique geodesic (up to rcparametrization) that connects P and Q. We can fix 
the affine parameter by requiring 2^(0) = P, x M (l) = Q. We now expand x^(X) in a power series in A around A = 0, 



xfiX n . 



(A2) 



The requirement that x^ (0) = P constrains to be equal to the coordinates of P in the chosen, arbitrary coordinate 
system. Given a set of three orthonormal space-like unit vectors (ei) p at P which are orthogonal to (eo)p, we can 
further write 



dx» 
~dX 



x=o 



■ x p (ei) 



p ■ 



(A3) 



Since (e;)^ are parallel transported along the central geodesic, they are uniquely defined along the geodesic once they 
are specified at one point; i.e., they are unique up to three Euler angles. Eq. (|A3[) defines the spatial Fermi coordinate 
x J F (recall that we have defined A through ^(A = 1) = Q, so that 5ijX F x F provides a measure for the spatial distance 
between points Q and P). 

In order to obtain the metric in the Fermi coordinate up to including 0{x F ), we also need quadratic (a^ ) and cubic 
(ag) order coordinate transformation. For that, we use the geodesic equation: 



d 2 x^ 
~dX 2 ~ 



^ dx a dx 13 
ap ~d~X~dX 



0. 



(A4) 



Rearranging the equation immediately yields 

d 2 x^ 



dX 2 



2a£ 



A=0 



dx a dx? 



1, 



ap dX dX 



x=o 



r/* — - -r M 

a 2 — 2 a/3 



_1 a/3 



( e i)p( e j)p x F x: F 



( e i)pi e j) P X F X F 



(A5) 



Applying one more derivative with respect to A to the geodesic equation yields through 



1 d 3 x>* 



<y>. 



6 dX 3 
1 

~ 6 



A=0 



1 d 
6dX 



^ dx a dx? 
af} 1X~dX 



A=0 



1 Q0,7 Zi (70L fi~f 



( e i)p( e j) p{&k)p x F x ?p x F- 



(A6) 
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Combining all, we find the coordinate transformation from FNC to general coordinates on a fixed x F hypersurfacc, 
up to third order in x l F : 



t p Z r 



l r 



p/A _ np/J- per 
1 a/3, 7 Z1 cra L £7 



(e l ) P (e j ) P (e k )lx i F x' F x k F . (A7) 



Here we have made explicit that (e v ) tl and T^p are always evaluated along the central geodesic at P. In App. lA l| we 
use this to derive the quadratic order FNC metric g^ v for a given global metric g^ v [Eq. (|A23|i ] . This result is also 
given in Eq. (66) of [l3| and Eq. (13.73) in [H]. In App. I A 2l we then apply this procedure to the perturbed FRW 
metric in synchronous gauge [Eq. ([!])] to obtain the Fermi coordinates and corrections oc xp to the metric. Readers 



familiar with Fermi coordinates and Eq. (|A23|) may want to skip to App. IA 2 



1. Metric in Fermi coordinate 

In this section we derive the metric in the Fermi coordinate by using the coordinate transformation we have found 
in Eq. (|A7[) . The final result will be that derived by [l~3|, Eq. (|A23[) . Under a coordinate transformation, the metric 
tensor transforms as 



F dx a dx 13 . . 

9uAxf) = 5-7; a „ g a p(x). 



dx F dx F 



(A8) 



Given the coordinate transform Eq. (IA7|) . we derive the partial derivatives to second order in xp, yielding 

( e i) a { e jY XpXp 



dx^ _ dx^(P) d i Id 

~^ + dx% [ei) XF ~2dxl 



dx F dx F 



dxf_ 
dxL 



1 a/3 



(eo)"-r^ ( ei ) a (e ) 



1 



pM 9T^ ~v a 

1 a/3,7 AL ap L -ya 



(eoy^ne^xpxp 



(A9) 



pM 9PM V a 

1 a/3,7 AL cj"y a/3 



pA 1 I opA 1 _ opA 1 p"" /LPA 1 P" 7 

1 a/3,7 "r" Z1 7 a,/3 zl cy af) ^ L ap L -fa 



(e i ) a (e j )( , (e l yx)?x' F + 2(e,) tt (e/(e t )V f 4 

( ei ) a (e,y(ei) 7 44. ( Al 0) 



where all unit vectors are evaluated at P. Note that dx t *(P)/dx F — (eo)p, since by definition (eo) M is the tangent 
vector to the central geodesic at P. Further, we have used that (ej) p and only depend on x° F , and that by 

construction, the basis vectors (ei) M are parallel transported along the central geodesic. This implies 







D 

DxJ 



dx F 



In Eq. (|A9|) we have also used 



(e0° = 



d 







dx° F a& 



_ pfi 



a/3,7 



(eo) 



(All) 



(A12) 



Finally, we need to take into account that the metric on the right-hand side of Eq. (|A8j) is evaluated at a point Q 
(specified by x F ) away from the central geodesic. We perform a Taylor expansion of g^ LU around P, 



+ 9aP,fi 



<5x M + -g a 0,ij,v 



6x»6x u + 0(Sx 3 ) 



(A13) 



where, from Eq. (IA7|) . 



8x^ — (ei)pX F n^ap i e i)p i e j)p x F x F 



(A14) 
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That is, up to second order in x % F , the metric at Q is given by 



g a p{Q) = g a p 



+ 9<*p„ 



(e.i) P x F + — \g a p,nv — gap.aF^] p (ei)p(ej) P x F x F . 



(A15) 



From now on, every instance of g pi/ , T°p, and (e„)' 1 will be evaluated at P, and we will omit P for brevity hereafter. 
Inserting Eq. (|A9[) . Eq. (|A10[) and Eq. (|A15[) into Eq. (|A8j) . and expanding to second order in x F yields the desired 
metric in FNC. At linear order, Eq. (|A8[) becomes 



F 

g^ 



(e M ) a -r« ( ei )"(e M r4 



p ^ 

F 



p-rMei)"WV F 

= (e M ) a (e„) g a p + (g a p. P - gap^ a - g a ^ a p p) (e M ) Q (e v ) p (e 8 ) p x 
= (e^)" {e v y g a p + g a p- p (e^) a {e v J (e l ) p x l F = t]^, 

where the last equality follows from the definition of the orthonormal tetrad at P, 

g^u{e a T{ep) v = r) a p, 

and the Levi-Civita connection, g^u-p = 0. 

Next, we calculate the quadratic correction to the metric in FNC, 



g a p + g a t3^(ek) K x F ] 



(A16) 
(A17) 



1 



r.g^.v.ap r.g{iv,a^ aj? ^g^,a^p v g^T a p v + 2g ll ~ j T a pY ua -\- g^aT^^Y ^ 



(e ) AI (e ) l/ (e i ) Q (e m ) /3 4^ 



7oi 



~^9iiv,aP ~ ^9fiiy,cr^ a fj ~ 9fij,a L fi v ~ 9iy-f,a^ ~ T^^W 1 a/3, fx ' 9v^ a a l .^cr 1 q^ 1 ^ 



-^gp,u,a/3 — -zQ 'fw ,oT \p ~ ,oXl [3 Stc^^F^ 



(e )' I (e 4 ) ly (eO Q (e m )^^F 



( ei )^(e,r(e ; ) Q (e m )^4^. 



Finally, using 



we have 



g^v.a g^v,a g<yv^ ' n.a .9^l '-^^ 



— "P CT _L P CT _1_ P ' P^ _L P ' P CT _L I P ' P CT _!_ P I P 

g^i/,ap — .9/jct 1 iyct.,3 i Scti/ 1 /ja,,3 i AW 1 crfi va ' ff?" 1 a/3 /JQ "+" 97cr I 1 ^^j 1 va < L vp L 



F 

(A18) 



(A19) 



(A20) 



(A21) 



(A22) 



By using Eqs. (|A21I) - (|A22[) . we write the partial derivatives of the metric as a function of the metric itself and the 
Christoffcl symbols, to obtain the final expression for the FNC metric at quadratic order: 

(A23) 



Sgoo = 


E)F I m 
JX OlOm x F J 'F 


Sgai = 


2 R F T l T m 
g rt O!im x F x F 


f«5 = 


_ I R F r l T m 

2 ) JX iljm d ' F x F 



Here, we have defined R^p^s as * ne Riemann tensor in FNC, 

R F aM s = (e a r(epr(e 7 r(e s ) x R^ K x, 
where the Riemann tensor is defined following the convention of (29| 

pM pM pM I pM ~pcr p/i per 

n apt ~ 1 cry,P L aP.j ^ L ap L a-y L <T~/ L ap- 



(A24) 
(A25) 
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2. Application to synchronous gauge metric 

We write the perturbed FRW metric Eq. ([I]) in terms of proper time t instead of conformal time rj, 

ds 2 = -dt 2 + a 2 (t)[d i:j + !>,, d.r'd.r 1 . (A26) 

For this metric, coordinate time coincides with proper time for a comoving observer (x l = const). Thus, without 
loss of generality we choose {(t, 0,0, 0)}t as central geodesic. Further, since the FNC time coordinate tp is given by 
the proper time along the central geodesic, we have t_F = t along the central geodesic (while £f 7^ t for x F 7^ 0). 
Correspondingly, the unit time vector is given by (eo)^ = (1, 0, 0, 0), and orthogonal spatial basis vectors are given by 
(see also (3l| ) 



(e k Y 



0, 



1 



5%k — 2^ ik 



(A27) 



The inverse metric is 



This leads to the following Christoffcl symbols, 



1 00 — 1 Qi — 1 00 — u 



1 



where here and throughout, dots denote derivatives with respect to t. The Riemann tensor is given by 



n 00 

pn 
11 Oi 



= [H + H 2 ] S im + -him + Hh 



2 \ hnrrL,i hni, 

& H [SnjSim Snm^ij] ~\~ ^ {^mn.ij ~\~ hij,nm him,nj ^jn.im) 
~l~ ^ ^ (foim^nj fo'ijfinm) ~t~ 2" ^^nj'^im ~t~ him&nj h%jS n m hnm^ij 



Wc then have 



RiOOm — gijR 3 



00m 



a 2 H + H 2 )S. 



lv 



H + H ] h 



(A28) 

(A29) 
(A30) 
(A31) 
(A32) 

(A33) 
(A34) 

(A35) 

(A36) 
(A37) 



Rlijm — ^ ^ [^ij^im ^Zm^ij] ~t~ c. (^m/,ij ~l~ hij,lm him,lj h>jl,im*) 



+ a 4 H 2 {h im 8ij + h[j5 im — hijSi m — h lm Sij) H — ( 



hljSi m + hi m S[j — hijSim — hl m 5ij 



Finally, the Riemann tensor in terms of FNC is given by 



1? 



!00m = ( e ( e o) ( e o) (e m ) -R^kA = (H + H 8i m + 



-hi m + Hhim 



R F 



iOim 



2a 

Rlijm = [dlj^im ~ ^Im^ij] + 2~^2 (^ m Mj ^ij,lm ~ hitnjj ~ ^jl.im) + ~ ( ^Ij^im + h{ m Slj — h>ijSl m — hl m S. 



1 



(A38) 

(A39) 
(A40) 

»)■ 

(A41) 
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Combining all with Eq. (|A23|) . we find that the metric in FNC is 

1 



F 

9oo 



F 

9oi 



F 



- 1 
1 



H + H 2 )rl 



I II I 



Hh 



im 



I pX p . 



— ^Vi/ljm Vm^ii ) XpXp 

H 2 

CpXp T p&ij 



(VjVj/imz + V/V m ft,y - ViV 3 h im - Vi\/ m hji) x l F x F 



H 



g | h[jX pX p ~h h{ rn x p x p h^j v p h[ 7n x pX p Sij 



(A42) 



Here, we define rp = SijX F x F and denote the partial derivative with respect to the FNC by V, = d/Ox F . Note that 
in Eq. (|A42|) . the derivative terms are already order x F , hence we can use Vj = {l/a)d/dx l here. We reiterate that 
Eq. (|A42|) is valid for any spatial metric perturbation hij , and thus also encompasses scalar cosmological perturbations 
written in synchronous-comoving gauge. 

It is also useful to have an explicit expression for the transformation from global coordinates x^ to Fermi coordinates 



Evaluating Eq. (1A7[) for the metric Eq. (|A26|) yields 







Sij 


2 kij 








~\hij 



^\k~^ x F x F + 0{x F ) 



(A43) 
(A44) 



This result is used in § VIIB of [1]. 

Finally, given the FNC metric Eq. (|A42|) , we can derive the motion of non-relativistic bodies (of momentum p 1 and 
mass to), which is governed by 



1 dp 1 
m dt 



F 
' 9oo,, 



(H 2 + H)x l F - V 4 * J 



1 fl 



9" = - - -h 



2 V 2 



Hh^j J x pX p 



(A45) 
(A46) 



This is the usual quasi-Newtonian equation of motion of a particle in an expanding Universe with peculiar potential 
^ F . The effective potential induces a tidal tensor given by 



tt 



0,0, - 

2 (-^hhn + Hhir 



OiOj 



SijW I x F x^ 



-him + Hhi m j - -StjTr ( -h hn + Hh h , 



If is traceless, the last term vanishes and we obtain 



tracclcss / 1 V tt'u 
tij = — [ -film + rlhlm 



(A47) 



(A48) 



Appendix B: Derivation of shear 
1. Shear statistics from tensor modes 



We express the tensor metric perturbation /i^ (x, 77) as 

MM)= J -^[4j^)h + (Kv)+efj(k)h^(k,r,) 
d 3 k 



p=-i,i 



(Bl) 
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where we have defined the helicity±2 polarization tensors and Fourier amplitudes through 



h±i = -(h+ T ih> 



(B2) 



In standard spherical coordinates, we have 



I cos 9 cos </> ± i sin ( 

m± = —j= (eg =F ie^) = —= cos 9 sin <ft =F i cos ( 
V2 V2 \ _ gin( 



V2 



(B3) 



where eg and are, respectively, the unit vectors of the polar and azimuthal angles. In order to make progress, we 
begin by evaluating the contribution of a single plane wave, assuming that k = kz. We have 



1 



4(k,n) = 4(k)m* T (nK T (n) = -(1 TPM )V 2 ^ 
ejf(k) = e\Mrj = (1 - M 2 )e i2 ^ 



e?,(k)m* (n)tf 



/'- 



V2 



(B4) 



where p = ±1 and /i = cos#. We will also use the notation k± — mlfe = — sin9 k/y/2. Using Eq. (|15p and Eq. (|20|) . 
we then have for the contribution to the shear 

( ±27 )(k,n)= J2 I -\[h P {K^) + h p {^ny^el + \CiP cM H- 2 h- 2 ^ 
p=-l,i I 



p=-i,i 



( 1 - -C liOcr oif - 2 a- 2 {9? + aHO^j h p (li,fj)e ikA 



1 



-(lTPM) e 



^ 2 (1 - P 2 ) 2 + (l - 2|) 4(1 - f)b TP) + ^(IT ptf 
4 X V X/ 2 2 X 



xe i2p \(k,^-x)e 



zk-n x 



(B5) 

Here, fj = 770 — x is the conformal time at emission inferred from the observed redshift, and all tilded quantities are 
evaluated at the source redshift. The next step will be to derive the spherical harmonic coefficients of the shear. 
Clearly, all factors involve e ±42< ^, so that only spherical harmonic coefficients with m = ±2 will be non-zero (this is of 
course a consequence of the choice k = kz). 



2. Spin-raising and lowering 



Let us consider the case of 27, and restrict to one circular polarization p = +1 first: 



27(k,n,+l) 



/ii(k,77o)e 



/.;// 



1 - \c lPcM H^ 2 ~a- 2 [d 2 + aHda}) Zn(k, fj)( 



(1 - M ) V 2 ^ 



dx 



*-^fc 2 (l - M 2 ) 2 + (l - 2| ) f (1 - ^ - 1) + Ul M ) 2 
4 X V X/ 2 x 2x 



(B6) 
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where we have defined x = k\, x = k\- Since this is a spin+2 quantity, we apply the spin-lowering o per ator twice to 
obtain a scalar. For this, we use that for functions that satisfy 8$ s f = im s f (see App. A of [lj and [32| for details), 



3 2 2fM) 



9 + T ^) 2 [(l-^ 2 ) 2 /(M, 



dfjb 



(B7) 



With m = 2, this yields 



e 2 27 (k,n,+l) 



1 / d 



4 V dn ' 1 - fi 2 

(l-/i 2 )(l- M ) 2 
2 



dfj, 1 — fi 2 



^(k,r? )e^ 

2 



a;=0 



1 - ^C lPcr0 H„ 2 a- 2 [dl + aHd^ hiQt, fj)e^ J, 



4 XX 



2 XX 2x 



1 / d 



4\ 3/i 1 - [i 2 
d X ( d 2 



^(k,r? )e^ 
2 (1-^ 2 ) 



z=0 
1 



1 — gCiPcrO 



F - 2 a- 2 {5? + off ft}}) /ii(k,r/)e^ }, 



i2tj> 



-x 2 {l-^ 2 ) 2 + -zx{l- P 2 )^-l) 
I (- l -x 2 {l M 2 ) 2 - ix{l - M 2 )(m - 1) + i(l - M) 2 



x e 42 ^i(k,7 7o -x)e MA1 



/'1 



(k,r? ) K?i(z)< 



2'=0 



+ 1 - oClP 



crO 



ff- 2 a- 2 {a 2 + affa rl } J hx(k, v)Qi(£) 



x dx 
x L 



Q 2 (x) + ±Q 3 (x) 
x 



^(k^o-x)^^ (l-^ 2 )e 420 , 



(B8) 



where we have turned powers of \x into powers of — id x and defined derivative operators Qi{x) as 



Qi(x) = 12- x 2 + 8xd x + x 2 d 2 x -i(8x + 2x 2 d x ) 

Q 2 {x) = - - [I4x 2 + x 4 + (40a; + 14a; 3 )<9 x + (50x 2 + 2x 4 )5 2 + 14x 3 d% + x A df\ - \i [4a: + x 3 + 6x 2 d x + x 3 d 2 ] 



h(x) = - [24 + 24a; 2 + x 4 + (96.x + 16x 3 )d x + {72x 2 + 2x i )d 2 + I6x 3 d 3 + x i d i x ] 



(B9) 



Although these operators are complicated, they will facilitate the connection with the convergence and rotation below. 
Note that, as expected, the real parts of Qi(x) only involve even powers of x (counting derivatives as well), while the 
imaginary parts involve odd powers only (where I111Q3 = 0). We now turn to p = — 1: 



27(k, n, -1) 



/i_i(k,ry )e 



X=0 



1 - ^C lPcr0 H^ 2 h- 2 [d 2 + ~aHd n }j /i-i(M) 



x-xx.2. 



4 X 



1 



^fc 2 (l - fj, ) + 1 - 2^ — (1 - M 2 )(- M - 1) + — (1 + fi) 



Xj 2x 



2x 



The appropriate spin-lowering operator is now for m = —2, i.e. 

(-8, 2/(1 - M 2 )) 2 = (8, + 2/(1 - M 2 )) 2 = + 2/(1 - m 2 )) 2 . 



(1 + ^) 2 e- i2 * 

e-^-iCk.fjo-xJe* 6 ". 

(BIO) 

(Bll) 
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Thus, 27(k, n, —1) is equal to 27(k, n, +1) [Eq. (|B6[) ] when changing /j, —> — /i, x —> —x in addition to hi —> h-i, 
—<j). Since Qi(—x) = Q*(x), we obtain 



8 2 27 (k,n,-l) 



o X 



+ ±Qt(?) 



(B12) 



Similarly, we can derive the corresponding expressions for -27, 



_ 27 (k,n,+l) 



(1 + /1) V 2< ^ 



X- XX, 2/1 2n2 



4 X 



X\ 



^ 2 (1 - n ) + 1 - 2± —(1 - M 2 )(- M - 1) + — (1 + fi) 



X) 2x 



2x 



e^iCk.fjo-xJe* 6 ", 
(B13) 



and correspondingly for p = —1, by acting twice with the spin-raising operator for m = ±2, (— =F 2/(1 — /i 2 ))- This 
immediately leads to 



! _ 27 (k,n,+l) = <^ - 



! _27(k,n,-l) 



fci(k,Tjo) (Qt(x)e ix ») 
X ^\q* 2 (x) + ^QI(x) 

xY x 
ft_i(k,i7o) (Qi(x)e J: 



x ^x 
X 



? 2 (x) + ^Q 3 (x) 
X 



e i ^ 1 (k,r 7o -x)|(l-M 2 )e i ^ 

1 - \c M H^ 2 ~a- 2 [d 2 + aHd^Yj fc_i(M)Qi(5) 
^^(Mo-xmi-M 2 )^ 20 - (B14) 



z=0 



This is in the desired form of Eq. (A17) in App. Al of 1.], with azimuthal harmonic index r = ±2. Since the shear 
±27 is a spin±2 quantity, and h±\ are two independent polarization states with power spectra Ph ±1 (k) = PTo(k)/8, 
we can apply Eq. (A24) in Q with s = 2, r = 2, N P = 2, and P h = P T o/8: 



fc 2 &P T0 (fc)|F^(fc)| 2 



T T (fc,77o) ReQi( 



x dx 
X 



ReQ 2 (a;) + ^ReQ 3 (:c) 
X 



ie=0 



1 — gCi/9 cr o-ff 5 



T T {k,r]Q - x) 



2~-2 



{9| + a#^}jT T (fc,r?)ReQi(x)^| 



(B15) 



X ^X 
X 



X 

Im<5 2 (2;) + — lmQ 3 (x) 



i^-T T {k,ri - x). 
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Again, fj = r/o — Xi 5 = and x = k\, x = k\. This completes the derivation of the angular power spectrum of 

E- and B-modes of the shear. The operators Qi when applied to spherical Bcssel functions can be simplified to yield 



ReQi(x) 
ImQi(x) 



RcQ 2 (x) J -^- 



ImQ 2 {x) — — 





1 


X 2 


X 2 




-i[ 


X 2 


X 



[(2a; 2 - I 2 - 31 - 2)jj(s) + 2xj;+i(a 
[2(Z + 2) - 7,(a:)-2a;i J+ i(i)] 



(l - 1) 



x 

;4 c;2 



(/4 _ 5 f + 4) M|) + 2 (Z 2 + I - 2 )^±iM 





+ w 


-1) 


(i- 


(I - 




2) jj 


(z) 



(i + l)(l-2) 



a 1 + 2)1 Ji(ar) 

RcQ3(a; )__ 

ImQg^^M =0 . 



(B16) 



In the limit of . 



for / = 2, we have 



ReQi(a:)-^- = - 
ReQ 3 {x) — — = -, 

X D 



(B17) 



while all other operators vanish in this limit for I = 2, and all operators vanish in this limit for I > 2. With this, we 
can easily verify that modes with k — > do not contribute to the quadrupole, as desired. As we let k — > 0, and thus 
x —> 0, we trivially have F^ B (k) — > 0, and 



^w-iHiJ+jf. t!UJ = -5 + 5=°- (B18) 

where we have used that Tx(k — > 0, rj) — > 1 (of course, we only need the fact that Tx(k — > 0, 77) — > const). 



Appendix C: Connection to convergence and rotation 



In this section, we cross-check Eq. (|B15I) with the angular power spectra of coordinate convergence and rotation, 
through the relations Eq. (|23|) . 



1. Angular power spectrum of coordinate convergence 

We begin with the general expression for the coordinate convergence k = -l/2dxAx\ derived in Q, restricted to 
synchronous-comoving gauge, and assuming a tranverse-traceless metric perturbation hij : 



3 1 f x 



-^ii-^u + (x-x)|vi|~h| 



3 f x 
4(^11)0 + J d X 



——d\\h\\ — —/in 

2 11 11 2x 11 



(CI) 
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This is equivalent to the expression used in [?J, but in a form more convenient for the comparison with the shear. 
Considering a single plane-wave perturbation oriented along the z-axis, with +1 circular polarization, we have 



K(n,k,+1) = -(1-^V 2 V^ h 1 (k, Vo )+ d X 

4 x=o 



dx^Ci-ziV'^MM), 

XX 



■-—xd x — — — 
2x 2 X 



(C2) 



where we have used Eq. (|B4[) . and turned z/xA; into xj\d x . understanding that the derivative only acts on e lx ^. In 
order to derive the multipole moments of k, we now use the relation (see App. A in [![) 



(C3) 



which yields 

«L(k) = - 



'2Z + l(Z + 2)' , 
4, (Z-2)! (47r),jma 



-/n(k,7jb)-^5- 



dxhi(k,v) 



-xd x 



2x 2 X 4 



3 + l i(I + 1) *-xU(*) 



XX J a; 



(C4) 



We thus obtain for the angular power spectrum of k: 
1 



if(fe) 



2tt 



'(1 + 2)1 
(l-2)\ 



k 2 dkP T0 (k)\F?(k)\ 2 

3l{x) 



;T T (k,Tj )- 



x=0 



d\ 
X 



1 



3 1 



T T {k, % - X) ~^xd x -- + -1(1 + 1) 1 



2 4 



1 Z(Z + 1) 
(Z + 2)(i-l) 



-(| + 2)(l-l)|T T (fc,i»)i* Ia 



(Z + 2)(Z - 1) / -^^(fc^o-x)^-^- ^ +4^ + 1) ( 1 



3 1 



x\ \ 3i{x) 



XJ ) x 



, (C5) 



where for the observer term wc have used lim^^o 3i( x )/ x2 = 1/15 for I = 2, and for I > 2. We can now simplify the 
operator applied to the spherical Bessel function, using the recurrence relation j[ — l/x ji — ji+i, yielding 



3i( x ) 



p+i-^o+!)]4-— 

2 ar x 

3i , lji+i 



(Z-2)(Z + 1)^ + 

x z 2 x 



(C6) 



Wc thus have 
F?(k) - 



1(1 + 1) 
(I + 2)(i-l) 



:T T (k,r]o)Si2 



[ — T T (k,rj - x) 
J X 



(Z + 2)(Z-1) 



( ,. 2){I + 1 )4 + 2 M}.^|±|2W 

a; 2 x J 4%(/ — 2J! a;" 1 



(C7) 



For comparison, the corresponding filter function for the i?-mode of the shear is (without metric shear and IA 
contributions, as discussed in § IV Ap 



F? E (k)= --T T (k, Vo )5 l2 - I ^ 



x dx 
o X 



(l + 2)(l-l) 



a + i) G -2)^ + 2 ^iM}_^([±|!Mj) 

x 2 x J %4(/-2)! x 2 



T T (k,m - x), 
(C8) 
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where we have used the x — > limit for KcQiji/x 2 for the observer term. Hence, 

F? E (k) = ^ (t t(?+ii 1) jr "W- (C9) 

Our result for the shear (without metric shear and intrinsic alignment) thus recovers the full-sky relation between 
shear F-(gradicnt-)modes and coordinate convergence Eq. (|2"3l) . 



2. Angular power spectrum of rotation 



Using the definition of the rotation Eq. (|2"2"T) . and the expression for Ax 1 ^ [Eq. (36) in [T[] restricted to synchronous 
gauge, we obtain 



"=2 
1 

~ 2 



dx e ljk n l (d J ± h k m )n r - 



e uh k ' 3 



(CIO) 



Note that eijkd 3 L n k = e^e^'n = 0, and n J '9j_ j/(x, it) = 0, and that pulling the derivatives inside the integrand yields 
a factor of x/x- This result agrees with Eq. (4) in (this ui is also equivalent to Vgfi as defined above Eq. (19) of 
|||). We now calculate the angular power spectrum of ui. Assuming as above a single plane wave with k = fez, we 
have 



jk h k Jh l h m = -ik [-2h 1 h 2 h+ + ((n 1 ) 2 - (n 2 ) 2 ) h> 
= — ik(l — /j 2 ) [— sin 2(f>h + + cos 24>h x ] e 



= fc(l-Ai 2 ) [/lie 



2i<j> 



since = h*. Thus, 



k) = ±£d X k [h ie 2i * - h- ie - 2i *] (i 



/x 2 )e lfe ^. 



(Gil) 



(C12) 



Note the relative minus sign between the two polarization states which shows that u> is parity-odd. In analogy with 
the derivation for k (see also App. Al in [J) we use Eq. (|C3[) to obtain the angular power spectrum of the rotation 
(it only contains "_B-modes"): 



i 

2^ 



k 2 dkP T0 (k)\Ff(k)\' 



™ s 4, /(, + 2) 



/ — xT T (k,Tio - x) 
J X 



3i{x) 



In Dodclson ct al. 



2 V (i-2)l 

an additional "metric shear" term was added to F"(k): 



x=kx 



(C13) 
(C14) 



f, w -> f, w + f; 



Fr' 



.MS 



(fe) = - 



'(1 + 2) 



2 (l + 2)(l - 1) V (Z-2)! 



(i - 1)— Ji-i{x) 



T T (k,fj). 



(C15) 



Note that in our convention, there is an overall minus sign for both Ff and F[ 
spectrum of the shear without observer, IA and FNC terms, Eq. (|B15[) with 



u,MS 



F? B (k)= + 2)(J-1) f X ^xT T (k, Vo -x) 

2 Jo X 



' (/ + 2)(/-l) 

Ki + i) 



Next, consider the B-modc power 

kip) 
x 2 

(C16) 
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where we have used Eq. (|B16|) and Eq. (|C14|) . We thus recover the relation Eq. (l23l) between the -B-modes of the 
shear and the rotation (without "metric shear"). The contribution of the FNC term to the -B-modes of the shear is 
given by 



f FNC (fc) = - ±T r (M)Im&(5;)^p = -1t t (M) 



(I -l) 3 ^-3i-i{x) 



1(1 + 1) 



(C17) 



showing that the "metric shear" contribution to u derived in [8[ agrees with the contribution of the FNC term to the 
shear B-modes. I 
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